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DERIVATION OF EIGENVECTORS FOR SPATIAL PROCESSING 
IN MIMO COMMUNICATION SYSTEMS 

Claim of Priority under 35 U.S.C. §119 
[0001] The present Application for Patent claims priority to Provisional Application 

No. 60/432,760 entitled "Derivation of Eigenvectors for Spatial Processing in MIMO 
Communication Systems" filed December 11, 2002, and assigned to the assignee hereof 
and hereby expressly incorporated by reference herein. 

BACKGROUND 

I. Field 

[0002] The present invention relates generally to data communication, and more 

specifically to techniques for deriving eigenvectors based on steered reference and used 
for spatial processing in multiple-input multiple-output (MIMO) communication 
systems. 

IL Background 

[0003] A MIMO system employs multiple (N T ) transmit antennas and multiple (N R ) 

receive antennas for data transmission. A MIMO channel formed by the Nt transmit 
and N R receive antennas may be decomposed into N s independent or spatial channels, 
where N s <min{AT r , N R }. Each of the N s independent channels corresponds to a 

dimension. The MIMO system can provide improved performance (e.g., increased 
transmission capacity and/or greater reliability) if the additional dimensionalities 
created by the multiple transmit and receive antennas are effectively utilized. 
[0004] In a wireless communication system, data to be transmitted is typically 

processed (e.g., coded and modulated) and then upconverted onto a radio frequency 
(RF) carrier signal to generate an RF modulated signal that is more suitable for 
transmission over a wireless channel. For a wireless MIMO system, up to Nt RF 
modulated signals may be generated and transmitted simultaneously from the N T 
transmit antennas. The transmitted RF modulated signals may reach the Nr receive 
antennas via a number of propagation paths in the wireless channel. The characteristics 
of the propagation paths typically vary over time due to various factors such as, for 
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example, fading, multipath, and external interference. Consequently, the RF modulated 
signals may experience different channel conditions (e.g., different fading and multipath 
effects) and may be associated with different complex gains and signal-to-noise ratios 
(SNRs). 

[0005] To achieve high performance, it is often necessary to estimate the response of 

the wireless channel between the transmitter and the receiver. For a MIMO system, the 
channel response may be characterized by a channel response matrix H , which includes 
N T N R complex gain values for N T N R different transmit/receive antenna pairs (i.e., one 
complex gain for each of the N T transmit antennas and each of the N R receive antennas). 
Channel estimation is normally performed by transmitting a pilot (i.e., a reference 
signal) from the transmitter to the receiver. The pilot is typically generated based on 
known pilot symbols and processed in a known manner (i.e., known a priori by the 
receiver). The receiver can then estimate the channel gains as the ratio of the received 
pilot symbols over the known pilot symbols. 

[0006] The channel response estimate may be needed by the transmitter to perform 

spatial processing for data transmission. The channel response estimate may also be 
needed by the receiver to perform spatial processing (or matched filtering) on the 
received signals to recover the transmitted data. Spatial processing needs to be 
performed by the receiver and is typically also performed by the transmitter to utilize 
the N s independent channels of the MIMO channel. 

[0007] For a MIMO system, a relatively large amount of system resources may be 

needed to transmit the pilot from the N T transmit antennas such that a sufficiently 
accurate estimate of the channel response can be obtained by the receiver in the 
presence of noise and interference. Moreover, extensive computation is normally 
needed to process the channel gains to obtain eigenvectors needed for spatial 
processing. In particular, the receiver is typically required to process the channel gains 
to derive a first set of eigenvectors used for spatial processing for data reception on one 
link and may further be required to derive a second set of eigenvectors used for spatial 
processing for data transmission on the other link. The derivation of the eigenvectors 
and the spatial processing for data transmission and reception are described below. The 
second set of eigenvectors typically needs to be sent back to the transmitter for its use. 
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As can be seen, a large amount of resources may be needed to support spatial processing 
at the transmitter and receiver. 
[0008] There is therefore a need in the art for techniques to more efficiently derive 

eigenvectors used for spatial processing in MIMO systems. 

SUMMARY 

[0009] Techniques are provided herein for deriving eigenvectors based on steered 

reference and used for spatial processing for data reception and transmission. A steered 
reference is a pilot transmission on only one spatial channel or eigenmode of a MIMO 
channel for a given symbol period, which is achieved by performing spatial processing 
with a steering vector for that eigenmode, as described below. The steered reference is 
used by a receiver to derive estimates of both a diagonal matrix Z of singular values 
and a unitary matrix U of left eigenvectors of the channel response matrix H , without 
having to estimate the MIMO channel response or perform singular value 
decomposition of H . 

[0010] The estimates of 2 and U may be used for matched filtering of data 

transmission received via a first link (e.g., the uplink). For a time division duplex 
(TDD) system, which is characterized by downlink and uplink channel responses that 
are reciprocal of one another, the estimate of U may also be used for spatial processing 
of data transmission on a second link (e.g., the downlink). 

[0011] In another aspect, a matrix U with orthogonal columns is derived based on the 

estimates of 2 and U . The orthogonalization of the columns of U may be achieved by 
various techniques such as QR factorization, minimum square error computation, and 
polar decomposition, all of which are described below. An orthogonal matched filter 
matrix M may then be derived based on the matrix U and the estimate of S . The 

matrix M may be used for matched filtering for the first link, and the matrix U may be 
used for spatial processing for the second link. 
[0012] Various aspects and embodiments of the invention are described in further detail 

below. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[0013] The features, nature, and advantages of the present invention will become more 

apparent from the detailed description set forth below when taken in conjunction with 
the drawings in which like reference characters identify correspondingly throughout and 
wherein: 

[0014] FIG. 1 shows a flow diagram of a process for deriving an orthogonal matched 

filter matrix M based on a steered reference; 
[0015] FIG. 2 shows a wireless communication system; 

[0016] FIG. 3 shows a frame structure for a TDD MIMO-OFDM system; 

[0017] FIG. 4 shows transmission of steered reference and data on the downlink and 

uplink for an exemplary transmission scheme; 
[0018] FIG. 5 shows a block diagram of an access point and a user terminal; and 

[0019] FIG. 6 shows a block diagram of the spatial processing performed by the access 

point and user terminal for data transmission on the downlink and uplink. 

DETAILED DESCRIPTION 
[0020] The word "exemplary" is used herein to mean "serving as an example, instance, 

or illustration." Any embodiment or design described herein as "exemplary" is not 
necessarily to be construed as preferred or advantageous over other embodiments or 
designs. 

[0021] The techniques described herein for deriving eigenvectors may be used for 

various MEMO communication systems. For example, these techniques may be used for 
single-carrier MIMO systems as well as multi-carrier MIMO systems. For clarity, these 
techniques are described below for a single-carrier MIMO system. 
[0022] The model for a single-carrier MIMO system may be expressed as: 

r = Hx + n , . Eq(l) 

where x is a "transmit" vector with N T entries for the symbols sent from the N T 
transmit antennas (i.e., x = [jc, x 2 ... x Nt ] t ); 
r is a "receive" vector with N R entries for the symbols received via the N R 
receive antennas (i.e., r = [r, r 2 ... r N ] r ); 
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H is an (N R xN T ) channel response matrix; 

n is a vector of additive white Gaussian noise (AWGN); and 

" r " denotes the transpose. 
The noise vector n is assumed to have components with zero mean and a covariance 
matrix of A„ = a 2 l , where I is the identity matrix and a 1 is the noise variance. 
[0023] The channel response matrix H may be expressed as: 







K2 ■ 




H = 

















Eq(2) 



where entry h ij9 for i€ {1 ... N R } and ye {1 ... N T }, is the coupling (i.e., complex 

gain) between the j-th transmit antenna and the i-th receive antenna. For simplicity, the 
channel response is assumed to be flat across the entire system bandwidth, and the 
channel response for each transmit/receive antenna pair can be represented by a single 
complex value h t . . Also for simplicity, the following description assumes that 

N R > N T , the channel response matrix H has full rank, and N s = N T < N R . 

[0024] The channel response matrix H may be "diagonalized" to obtain the N T 

independent channels, which are also referred to as spatial channels or eigenmodes. 
This diagonalization may be achieved by performing either singular value 
decomposition of the channel response matrix H or eigenvalue decomposition of the 
correlation matrix of H, which is H"H, where " H " denotes the conjugate transpose. 
For clarity, singular value decomposition is used for the following description. 

[0025] The singular value decomposition of the channel response matrix H may be 

expressed as: 

H = U£V " , Eq(3) 
where U is an (N R xN R ) unitary matrix whose columns are left eigenvectors of H ; 

Z is an (N R xN T ) diagonal matrix of singular values of H, which is 

Z = diag (p u a 2 2 ... <t Nt Nt ) ; and 
V is an (Af r xN T ) unitary matrix whose columns are right eigenvectors of H . 
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A unitary matrix M is characterized by the property M W M = I , which means that the 
columns of the unitary matrix are orthogonal to one another and the rows of the matrix 
are also orthogonal to one another. The columns of the matrix V are also referred to as 
steering vectors. Singular value decomposition is described in further detail by Gilbert 
Strang in a book entitled "Linear Algebra and Its Applications," Second Edition, 
Academic Press, 1980. 

[0026] Spatial processing may be performed by both the transmitter and the receiver to 

transmit data on the N T spatial channels of the MEMO channel. The spatial processing 
at the transmitter may be expressed as: 

x = Vs , Eq(4) , 

where s is a "data" vector with up to N T non-zero entries for data symbols to be 
transmitted on the Nt spatial channels. The transmit vector x is further processed and 
then transmitted over the MIMO channel to the receiver. 

[0027] The received transmission at the receiver may be expressed as: 

r = Hx + n = HVs + n = U£V"Vs + n = UZs + n , Eq (5) 

where all the terms are defined above. 
[0028] The spatial processing at the receiver to recover the data vector s may be 

expressed as: 

s = GMr = GZ r U"r = GS r U" (UZs + n) = s + n , Eq (6) 

where s is the data vector; 

s is an estimate of the data vector s ; 

M is an (N T xN R ) matched filter matrix, which is M = X r U" ; 
G is an (N T xN T ) scaling matrix, which is 
G = diag (1 / al 1 / o\ 2 ... U<t 2 Nt Nt ) ; and 

n is the post-processed noise, which is n= GZ r U w n . 
The spatial processing by the receiver is often referred to as matched filtering. Since 
M = 1?}J H and since the columns of U are left eigenvectors of H , the columns of M r 
are conjugated left eigenvectors of H scaled by the singular values in Z . 
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[0029] As shown in equation (6), the receiver needs good estimates of the matrices S 

and U in order to perform the matched filtering to recover the data vector s. The 
matrices X and U may be obtained by transmitting a pilot from the transmitter to the 
receiver. The receiver can then estimate the channel response matrix H based on the 
received pilot and perform the singular value decomposition of this estimate, as shown 
in equation (3), to obtain the matrices 2 and U . However, as noted above, a large 
amount of resources may be needed to transmit this pilot and to perform the singular 
value decomposition. 
I. Steered Reference 

[0030] In an aspect, a steered reference is transmitted by the transmitter and used by the 

receiver to derive estimates of the matrices S and U , which are needed for matched 
filtering. The steered reference is a pilot transmission on only one spatial channel or 
eigenmode for a given symbol period, which is achieved by performing spatial 
processing with a steering vector for that eigenmode. The receiver can then estimate the 
matrices Z and U based on the steered reference, without having to estimate the 
MIMO channel response or perform the singular value decomposition. 

[0031] A steered reference sent by the transmitter may be expressed as: 

*sr,m = v OT • p , for m e {1 ... N T } , Eq (7) 

where x sr m is the transmit vector for the steered reference for the m-th eigenmode; 
v m is the right eigenvector of H for the m-th eigenmode; and 
p is a pilot symbol transmitted for the steered reference. 
The eigenvector v m is the m-th column of the matrix V , where V = [y, y 2 ... \_ Nt ] . 

[0032] The received steered reference at the receiver may be expressed as: 

r sr , m = Hx sr , m + n = Hy m/ ? + n = UZV" v mj p + n = u m (7 m p -+ n , Eq (8) 
where r sr m is the receive vector for the steered reference for the m-th eigenmode; and 
a m is the singular value for the m-th eigenmode. 

[0033] As shown in equation (8), at the receiver, the received steered reference in the 

absence of noise is equal to u m cr m p , which is the known pilot symbol p transformed by 

u m cr m . The eigenvector u m is the m-th column of the matrix U , and the singular value 
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a m is the m-th diagonal element of the matrix Z . The receiver can thus obtain an 
estimate of u m c m based on the steered reference sent by the transmitter. 
[0034] Various techniques may be used to process the received steered reference to 

obtain estimates of u m and a m . In one embodiment, to obtain an estimate of , the 

receive vector for the steered reference sent on the m-th eigenmode is first 

multiplied with the complex conjugate of the pilot symbol, p* . The result may then be 
integrated over multiple steered reference symbols received for each eigenmode m to 
obtain the estimate of u m a m . A row vector m m may be defined to be equal to the 

conjugate transpose of the estimate of u m cr m (i.e., m m = & m u" ). Each of the Nr entries 
of the vector m m is obtained based on a corresponding one of the Nr entries of the 
vector . 

[0035] The row vector m m for the m-th eigenmode includes estimates of both n m and 

a m , and may thus be referred to as a scaled vector. Since eigenvectors have unit power, 
the singular value a m may be estimated based on the received power of the steered 
reference, which can be measured for each eigenmode. In particular, the singular value 
estimate & m may be set equal to the square root of the power for the vector r sr m , 

divided by the magnitude of the pilot symbol p. The vector m m may be scaled by l/<7 m 
to obtain the eigenvector u m . 
[0036] In another embodiment, a minimum mean square error (MMSE) technique is 

used to obtain an estimate of u m based on the receive vector r srm for the steered 

reference. Since the pilot symbol p is known, the receiver can derive an estimate of u m 
such that the mean square error between the recovered pilot symbol p (which is 
obtained after performing matched filtering on the receive vector r STm ) and the 

transmitted pilot symbol p is minimized. 
[0037] The steered reference is transmitted for one eigenmode at a time (i.e., one 

eigenmode for each symbol period of steered reference transmission). The steered 
reference for all Nt eigenmodes may be transmitted in various manners. In one 
embodiment, the steered reference is transmitted for one eigenmode for each frame, 
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where a frame is an interval of data transmission for the system and is defined to be of a 
particular time duration (e.g., 2 msec). For this embodiment, the steered reference for 
multiple eigenmodes may be transmitted in multiple frames. In another embodiment, 
the steered reference is transmitted for multiple eigenmodes within one frame. This 
may be achieved by cycling through the Nt eigenmodes in Nt symbol periods. For both 
embodiments, the n-th steered reference symbol may be expressed as: 

x S r,m(") = v (/imod ^ ]+1 -p , for ne{l...L}, Eq (9) 

where n is an index for either symbol period or frame number and L is the number of 
steered reference symbols to be transmitted. Multiple steered reference symbols may be 
transmitted for each eigenmode m to allow the receiver to obtain more accurate estimate 

[0038] The receiver is able to obtain the row vector m m for each of the Nt eigenmodes 

based on the received steered reference for that eigenmode. The row vectors for all 
Nt eigenmodes may be used to form an initial matched filter matrix M, where 

M = [mi lib — Oi/v r ] r anc * M = f^u" • The matrix M may be used for matched 
filtering by the receiver, as shown in equation (6), to recover the transmitted data vector 
s . 

[0039] The steered reference is sent for one eigenmode at a time and may be used by 

the receiver to obtain the matched filter vector m m for that eigenmode. Since the N T 

matched filter vectors m m of the matrix M are obtained individually and over different 
symbol periods, and due to noise and other sources of degradation in the wireless 
channel, the Nt vectors m m of the matrix M are not likely to be orthogonal to one 
another. If the N T vectors m m are thereafter used for matched filtering of a received 
data transmission, then any errors in orthogonality among these vectors would result in 
cross-talk between the individual symbol streams sent on the N T eigenmodes. The 
cross-talk may degrade performance. 
II. Eigenvector Orthogonalization 
[0040] In another aspect, to improve performance, an enhanced matched filter matrix 

M is derived based on the steered reference and has row vectors that are forced to be 
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orthogonal to one other. The orthogonalization of the row vectors of M may be 
achieved by various techniques such as QR factorization, minimum square error 
computation, and polar decomposition. All of these orthogonalization techniques are 
described in detail below. Other orthogonalization techniques may also be used and are 
within the scope of the invention. 
1. OR Factorization 

/v r 

[0041] QR factorization decomposes the transposed initial matched filter matrix, M , 

into an orthogonal matrix Q p and an upper triangle matrix R F . The matrix Q p forms 

/v T 

an orthogonal basis for the columns of the matrix M (i.e., the rows of M), and the 
diagonal elements of the matrix R F give the length of the components of the columns of 

M in the directions of the respective columns of Q p . The matrices Q p and R F may 

be used to derive an enhanced matched filter matrix M F . 
[0042] The QR factorization may be performed by various methods, including a Gram- 

Schmidt procedure, a householder transformation, and so on. The Gram-Schmidt 
procedure is recursive and may be numerically unstable. Various variants of the Gram- 
Schmidt procedure have been devised and are known in the art. The "classical" Gram- 

Schmidt procedure for orthogonalizing the matrix M is described below. 

[0043] For QR factorization, the matrix M may be expressed as: 

M r =Q F R F , Eq(10) 

where Q p is an (N R xN R ) orthogonal matrix; and 

R F is an (N R xN r ) upper triangle matrix with zeros below the diagonal and 
possible non-zero values along and above the diagonal. 
[0044] The Gram-Schmidt procedure generates the matrices Q p and R F column-by- 

column. The following notations are used for the description below: 
= [Q Q • * Q 1 > where q . is the y-th column of Q c ; 

q i j is the entry in the i-th row and j-th column of Q p ; 

Q F =[Q, Q 2 — ], where q^. is the y-th column of Q p ; 
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r. j is the entry in the z-th row and 7-th column of R F ; 

~ T ~ T 

M - [m, m 2 ... ] , where m y is the 7-th column of M ; and 

~ T 

m i j is the entry in the i-th row and 7-th column of M . 
[0045] The first column of Q F and R F may be obtained as: 



Li = II Oli II = 



Zha I 2 



1/2 



,and Eq(ll) 



q =— III, ■ Eq(12) 

The first column of R F includes one non-zero value r hl for the first row and zeros 
elsewhere, where r x x is the 2-norm of m, . The first column of Q p is a normalized 

- T 

version of the first column of M , where the normalization is achieved by scaling each 
entry of ifi! with the inverse of r u . 

[0046] Each of the remaining columns of Q p and R F may be obtained as follows: 

FOR 7 = 2, 3 ... N T 
FOR i = 2 ... 7-1 

r U = %^J Eq(13) 
% =™>-Z^ 3, Eq(14) 



1=1 



rjj = \\%\\ Eq(15) 

q. =— q. . Eq(16) 

[0047] The Gram-Schmidt procedure generates one column at a time for the matrix Q p . 

Each new column of Q F is forced to be orthogonal to all prior-generated columns to the 

left of the new column. This is achieved by equations (14) and (16), where the 7-th 
column of Q F (or q^ ) is generated based on , which in turn is generated based on the 

j-th column of M (or m y ) and subtracting out any components in m y pointing in the 
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direction of the other (j -1) columns to the left of m ; . The diagonal elements of R F 
are computed as the 2-norm of the columns of Q F (where q j = m, ), as shown in 
equation (15). 

[0048] Improved performance may be attained by ordering the matrix M based on the 

singular value estimates before performing the QR factorization. The initial singular 
value estimates a m , for m e {1 ... N T } , for diagonal matrix S may be computed as the 

~ T 

2-norm of the columns of M , as described below. The initial singular value estimates 
may then be ordered such that {a x > a 2 > ... > a Nr } , where a x is the largest singular 

value estimate and a Nr is the smallest singular value estimate. When the initial 

singular value estimates for the diagonal matrix 2 are ordered, the columns of the 

- T - T 

matrix M are also ordered correspondingly. The first or left-most column of M 

would then be associated with the largest singular value estimate and the highest 

received SNR, and the last or right-most column of M would be associated with the 
smallest singular value estimate and the lowest received SNR. For the QR factorization, 

the initial singular value estimates may be obtained as the 2-norm of the columns of 

~ r ~ t 

M and used for ordering the columns of M . The final singular value estimates are 

obtained as the 2-norm of the columns of Q p , as described above. The steered 

reference may also be transmitted in order (e.g., from the largest eigenmode to the 
smallest eigenmode), so that the singular value estimates are effectively ordered by the 
transmitter. 

~ T 

[0049] If the columns of M are ordered based on decreasing values of their associated 

singular value estimates, then the columns/eigenvectors of Q p are forced to be 

orthogonal to the first column/eigenvector with the best received SNR. This ordering 
thus has the beneficial effect of rejecting certain noise components of each of the 
remaining eigenvectors of Q p . In particular, the y-th column of Q F (or q^) is 

generated based on the y-th column of M (or m ; ), and noise components in m ; that 
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point in the direction of the j - 1 eigenvectors to the left of q^ (which are associated 

with higher received SNRs) are subtracted from m y to obtain q^ . The ordering also 

has the beneficial effect of improving the estimates of eigenvectors associated with 
smaller singular values. The overall result is improved performance, especially if the 
orthogonalized eigenvectors of Q F are used for spatial processing for data transmission 

on the other link, as described below. 
[0050] The enhanced orthogonal matched filter M F obtained based on QR factorization 

may then be expressed as: 

M F =Q F R F , Eq(17) 
where R F includes only the diagonal elements of R F (i.e., the elements above the 
diagonal are set to zeros). The diagonal elements of R F and R F are estimates of the 
singular values of H. Since M = ? r U" and M F = RfQ f > the following substitutions 
may be made: R F = Z and Q p ~ U* , where denotes the complex conjugate. 
2. Mean Square Error Computation and Polar Decomposition 
[0051] The initial matched filter matrix M may also be orthogonalized based on a 

particular optimality criterion. One possible criterion is to minimize a measure of the 
squared error between the matrix M and an "optimum" matched filter with the desired 
orthogonality properties. This may be expressed as: 

minimize || M-S r Q p || F subjectto QpQ P =I, Eq (18) 

where || X \\ F is the Frobenius norm of X , and is given as: 




Eq(19) 



The condition Q F Q p =1 ensures that Q p is a unitary matrix, which would mean that 
the columns of Q p are orthogonal to one another and the rows of Q p are also 
orthogonal to one another. Equation (18) results in an optimum matched filter ? r Q P 
that is the best fit to the measured data given by the matrix M . 
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[0052] The solution to equation (18) can be obtained from the known solution to the 

orthogonal Procrustes problem. This problem asks the question - given two known 
matrices A and B, can a unitary matrix Q p be found that rotates B into A. The 

problem may be expressed as: 

minimize || A - BQ p || F subject to Q p Q p = I . Eq (20) 

[0053] The solution to the Procrustes problem can be obtained as follows. First, a 

matrix C P is defined as C P =B H A. The singular value decomposition of C p is then 
given as C P =U P ZpVp or U P C P V P =Z P . The unitary matrix Q p that solves the 
minimization problem shown in equation (20) is then given as: 

Q P = UpYp y . Eq(21) 

The derivation and proof for equation (21) is described by G. H. Golub and C. F. Van 
Loan in "Matrix Computation", Third Edition, Johns Hopkins University Press, 1996. 
[0054] The solution for equation (20), which is shown in equation (21), is related to the 

polar decomposition of the matrix C . This polar decomposition is given as: 

C P = Z P P P , Eq(22) 
where Z p is a unitary matrix, which is given as Z p = U p V p ; 

U P is a matrix of left eigenvectors of C p that spans the column space of C P 
(i.e., Up is equal to U P or a sub-matrix of U p depending on the 
dimension of C P ); 

P p is a Hermitian symmetric positive semi-definite matrix, which is given as 
Pp=Ypl P Vp ;and 

S p is a square matrix of singular values of C P with dimension equal to the 
number of columns of C p . 
[0055] Polar decomposition can thus be performed on the matrix C P to obtain the 

unitary matrix Z p , which may be equal to either Q p or a sub-matrix of Q p depending 

on the dimension of C P . It can be shown that the matrix Z P is the optimal result to the 
minimization problem shown in equation (20). 
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[0056] Algorithms for direct computation of polar decomposition are described by P. 

Zielinski and K. Zietak in "The Polar Decomposition — Properties, Applications and 

Algorithms," Annals of the Polish Mathematical Society, 38 (1995), and by A. A. 

Dubrulle in "An Optimum Iteration for the Matrix Polar Decomposition," Electronic 

Transactions on Numerical Analysis, Vol. 8, 1999, pp. 21-25. 
[0057] The solution for the optimum matched filter expressed in equation (18) may be 

obtained based on the solution to the orthogonal Procrustes problem described above. 

This may be achieved by equating M to A and £ to B. For the computation, an 
estimate of the singular values, X , may be obtained as the 2-norm of the columns of 

* T ~ 

M and used in place of S . The diagonal elements of S may be expressed as: 

il/2 



= II m, ll = 



for ie{l ... N T }. Eq(23) 



It can be shown that the use of Z in the computation for Q p results in nearly un- 

measurable degradation in performance relative to the use of the exact singular values in 
Z. 

[0058] A matrix C M may then be defined as: 

C M =|M . Eq(24) 
The singular value decomposition of the matrix C M is then given as: 

C M = U M Z M V" M or UmC m V m = E M . Eq (25) 

The unitary matrix Q M that solves the minimization problem shown in equation (18) is 
then given as: 

Q m =UmYm=u" . ^(26) 

An enhanced orthogonal matched filter M M , which is the solution to the minimization 
problem in equation (18), may then be expressed as: 

M M =l r u"=l T U M V£ . Eq(27) 
[0059] Alternatively, the polar decomposition of C M may be performed as described 

above, which may be expressed as: 

C M =Z M P M . Eq(28) 
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The unitary matrix Q M that solves the minimization problem shown in equation (18) 
may then be given as: 

Q M =Z M =U" . Eq(29) 
[0060] The enhanced orthogonal matched filter M M may then be expressed as: 

M M =l r Z M . Eq(30) 
It can be shown that the matrix Z M from the polar decomposition is the optimal result 
for the matrix Q M for the minimum square error computation (i.e., Q M = Z M ). Thus, 
the polar decomposition and minimum square error computation both yield the same 
orthogonal matched filter M M . 
[0061] FIG. 1 shows a flow diagram of an embodiment of a process 100 for deriving an 

orthogonal matched filter matrix M based on a steered reference. Initially, the receiver 
receives and processes the steered reference to obtain an estimate of u m a m for each of 
multiple eigenmodes of H (step 112). This processing may be performed as described 
above. An initial matched filter matrix M is then formed whose rows m m , for 
me{l ... N T ], are derived based on the estimates of u m a m . The orthogonal matched 

filter matrix M may then be obtained from the initial matched filter matrix M using 
any one of the orthogonalization techniques described above. 
[0062] For the QR factorization technique, the matrix M is factorized to obtain the 

matrices Q p and R F (step 122). The orthogonal matched filter matrix M is then 

obtained as shown in equation (17) (step 124) and the singular value estimates X are 
obtained as the diagonal elements of R F (step 126). 
[0063] For the minimum square error technique, estimates of the singular values, L , are 

obtained as the 2-norm of the columns of M (step 132). The matrix C M is then 
computed as shown in equation (24) (step 134). Singular value decomposition of C M is 
next computed as shown in equation (25) (step 136). The orthogonal matched filter 
matrix M is then obtained as shown in equation (27) (step 138). 
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[0064] For the polar decomposition technique, estimates of the singular values, X , are 

obtained as the 2-norm of the columns of M (step 142). The matrix C M is then 
computed as shown in equation (24) (step 144). Polar decomposition of C M is next 
computed as shown in equation (28) (step 146). The orthogonal matched filter matrix 
M is then obtained as shown in equation (30) (step 148). 

[0065] The orthogonal matched filter matrix M may thereafter be used to perform 

matched filtering of a received data transmission (step 150). 
[0066] The orthogonalization of the matched filter matrix provides several benefits. 

First, the use of an orthogonal matched filter matrix M avoids cross-talk between the 

eigenmodes of H. The derivation of the initial matched filter matrix M piecemeal 

based on the steered reference does not guarantee that the eigenvectors of M are 
orthogonal. The lack of orthogonality results in performance degradation. The 
orthogonalization of the matched filter matrix avoids this performance degradation. 

[0067] Second, QR factorization can improve the quality of the eigenvectors associated 

with smaller singular values. Without QR factorization, the quality of the estimates of 
the eigenvectors is not constant, and the estimates of the eigenvectors associated with 
smaller singular values are likely to be lower in quality. QR factorization can improve 
the quality of the eigenvectors associated with smaller singular values by rejecting 
certain noise components, as described above. Polar decomposition may have similar 
effect, but not in the direct way as QR factorization. 

[0068] Third, orthogonalization may reduce the amount of resources needed to transmit 

the steered reference. If orthogonalization is not performed, then high quality estimates 
of Z and U would be needed to ensure low cross-talk among the eigenmodes. A 
longer transmission period would then be needed for the steered reference for the 
eigenvectors associated with smaller singular values to ensure that the desired quality is 
obtained. High quality estimates of 2 and U would thus require a longer period of 
transmission for the steered reference (which would consume more valuable system 
resources) and a longer integration period for the steered reference at the receiver 
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(which may result in longer delay for data transmission). Orthogonalization can provide 
the desired performance without the need for high quality estimates of Z and U . 

III. MIMO-OFDM System 
[0069] The techniques for deriving eigenvectors used for spatial processing are now 

described for an exemplary wideband MIMO communication system that employs 
orthogonal frequency division multiplexing (OFDM). OFDM effectively partitions the 
overall system bandwidth into a number of (N F ) orthogonal subbands, which are also 
referred to as tones, frequency bins, or frequency subchannels. With OFDM, each 
subband is associated with a respective subcarrier upon which data may be modulated. 
For a MIMO-OFDM system, each subband may be associated with multiple 
eigenmodes, and each eigenmode of each subband may be viewed as an independent 
transmission channel. 

[0070] For OFDM, the data or pilot to be transmitted on each usable subband is first 

modulated (i.e., mapped to modulation symbols) using a particular modulation scheme. 
One modulation symbol may be transmitted on each usable subband in each symbol 
period. A signal value of zero may be sent for each unused subband. For each OFDM 
symbol period, the modulation symbols for the usable subbands and zero signal values 
for the unused subbands (i.e., the modulation symbols and zeros for all N F subbands) are 
transformed to the time domain using an inverse fast Fourier transform (IFFT) to obtain 
a transformed symbol that comprises N F time-domain samples. To combat inter-symbol 
interference (ISI) caused by frequency selective fading, a portion of each transformed 
symbol is often repeated (which is often referred to as adding a cyclic prefix) to form a 
corresponding OFDM symbol. The OFDM symbol is then processed and transmitted 
over the wireless channel. An OFDM symbol period, which is also referred to as a 
symbol period, corresponds to the duration of one OFDM symbol. 

[0071] For this exemplary system, the downlink and uplink share a single frequency 

band using time-division duplex (TDD). For a TDD MIMO-OFDM system, the 
downlink and uplink channel responses may be assumed to be reciprocal of one another. 
That is, if H(£) represents a channel response matrix from antenna array A to antenna 
array B for subband k, then a reciprocal channel implies that the coupling from array B 
to array A is given by H T (k) . 
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[0072] FIG. 2 shows a wireless communication system 200 that includes a number of 

access points (APs) 210 that communicate with a number of user terminals (UTs) 220. 
(For simplicity, only one access point is shown in FIG. 2.) An access point may also be 
referred to as a base station or some other terminology. Each user terminal may be a 
fixed or mobile terminal and may also be referred to as an access terminal, a mobile 
station, a remote station, a user equipment (UE), a wireless device, or some other 
terminology. Each user terminal may communicate with one or possibly multiple 
access points on the downlink and/or the uplink at any given moment. The downlink 
(i.e., forward link) refers to transmission from the access point to the user terminal, and 
the uplink (i.e., reverse link) refers to transmission from the user terminal to the access 
point. The channel response between each access point and each user terminal may be 
characterized by a set of channel response matrices H(fc), for ke K where K 
represents the set of all subbands of interest (e.g., the usable subbands). 

[0073] In the following description for a pair of communicating access point and user 

terminal, it is assumed that calibration has been performed to account for differences 
between the transmit and receive chains of the access point and the user terminal. The 

results of the calibration are diagonal matrices K ap (k) and K ut (k) , for k e K , to be 
used at the access point and the user terminal, respectively, on the transmit path. A 
"calibrated" downlink channel response, H cdn (&), observed by the user terminal and a 

"calibrated" uplink channel response, H cup (£) , observed by the access point may then 
be expressed as: 

H cdn (fc) = H dn (^)K ap (^) , forfcetf, Eq(31a) 
EU (*) = H up (*)K m (*) , for k G K , and Eq (3 lb) 

H^W-HLp^) , for ke K , Eq (31c) 

where H^^t) = R ut (k)H(k)T ap (k) is the "effective" downlink channel response, which 

includes the responses of the transmit chain T ap (fc) at the access point 

and the receive chain R m (fc) at the user terminal; 
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Hup(fc) = R ap (£)H r (£)T ut (fc) is the "effective" uplink channel response, which 
includes the responses of the transmit chain T ut (fc) at the user terminal 
and the receive chain R ap (fc) at the access point; and 
H(fc) is an (N^ xN ap ) channel response matrix between the N ap antennas at the 
access point and the N ut antennas at the user terminal. 
If calibration is not performed, then the matrices K ap (&) and K ul (£), for k e K , are 
each set to the identity matrix I . 
[0074] FIG. 3 shows an embodiment of a frame structure 300 that may be used for a 

TDD MIMO-OFDM system. Data transmission occurs in units of TDD frames, with 
each TDD frame covering a particular time duration (e.g., 2 msec). Each TDD frame is 
partitioned into a downlink phase and an uplink phase. The downlink phase is further 
partitioned into multiple segments for multiple downlink transport channels. In the 
embodiment shown in FIG. 3, the downlink transport channels include a broadcast 
channel (BCH), a forward control channel (FCCH), and a forward channel (FCH). 
Similarly, the uplink phase is partitioned into multiple segments for multiple uplink 
transport channels. In the embodiment shown in FIG. 3, the uplink transport channels 
include a reverse channel (RCH) and a random access channel (RACH). 
[0075] In the downlink phase, a BCH segment 310 is used to transmit one BCH 

protocol data unit (PDU) 312, which includes a beacon pilot 314, a MIMO pilot 316, 
and a BCH message 318. The beacon pilot is transmitted from all access point antennas 
and is used by the user terminals for timing and frequency acquisition. The MIMO pilot 
is transmitted from all access point antennas with different orthogonal codes and is used 
by the user terminals for channel estimation. The BCH message carries system 
parameters for the user terminals in the system. An FCCH segment 320 is used to 
transmit one FCCH PDU, which carries assignments for downlink and uplink resources 
and other signaling for the user terminals. An FCH segment 330 is used to transmit one 
or more FCH PDUs 332. Different types of FCH PDU may be defined. For example, 
an FCH PDU 332a includes only a data packet 336a, and an FCH PDU 332b includes a 
downlink steered reference 334b and a data packet 336b. 
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[0076] In the uplink phase, an RCH segment 340 is used to transmit one or more RCH 

PDUs 342 on the uplink. Different types of RCH PDU may also be defined. For 
example, an RCH PDU 342a includes an uplink steered reference 344a and a data 
packet 346a. An RACH segment 350 is used by the user terminals to gain access to the 
system and to send short messages on the uplink. An RACH PDU 352 may be sent 
within RACH segment 350 and includes an uplink steered reference 354 and a message 
356. 

[0077] For the embodiment shown in FIG. 3, the beacon and MIMO pilots are sent on 

the downlink in the BCH segment in each TDD frame. A steered reference may or may 
not be sent in any given FCH/RCH PDU. A steered reference may also be sent in an 
RACH PDU to allow the access point to estimate pertinent vectors during system 
access. 

[0078] For simplicity, the following description is for a communication between one 

access point and one user terminal. The MIMO pilot is transmitted by the access point 
and used by the user terminal to obtain an estimate of the calibrated downlink channel 

response, H cdn (£), for ke K . The calibrated uplink channel response may then be 

~ T 

estimated as H^fc) = H cdn (fc). Singular value decomposition may be performed to 
diagonalized the matrix H^CA:) for each subband, which may be expressed as: 

H cup (*) = U ap (*)2(*)Y* (*) , for k G K , Eq (32) 

where U ap (£) is an (N ap x N ap ) unitary matrix of left eigenvectors of H cup (£) ; 

£(£) is an (N^xN^) diagonal matrix of singular values of H^fc) ; and 

V ul (fc) is an (N^ xN u! ) unitary matrix of right eigenvectors of H^CA;) . 
[0079] Similarly, the singular value decomposition of the estimated calibrated downlink 

channel response matrix, H c dn (fc) , may be expressed as: 

Hcdn(^) ~ Yl (k)±(k)vl p (*) , for k G K , Eq (33) 

where the matrices V ut (&) and U ap (£) are unitary matrices of left and right 
eigenvectors, respectively, of H cdn (£) . 
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[0080] As shown in equations (32) and (33), the matrices of left and right eigenvectors 

for one link are the complex conjugate of the matrices of right and left eigenvectors, 
respectively, for the other link. For simplicity, reference to the matrices U ap (it) and 

V m (fc) in the following description may also refer to their various other forms (e.g., 

V ut (*) may refer to V ul (*), \ ut (k), V^ t (*), and v"(fc)). The matrices U ap (*) and 

V ut (£) may be used by the access point and user terminal, respectively, for spatial 

processing and are denoted as such by their subscripts. The matrix includes 
singular value estimates that represent the gains for the independent channels (or 
eigenmodes) of the channel response matrix H(&) for the Jt-th subband. 
[0081] The singular value decomposition may be performed independently for the 

channel response matrix H cup (fc) for each of the usable subbands to determine the Ns 
eigenmodes for the subband. The singular value estimates for each diagonal matrix 
Z(fc) may be ordered such that {a l (k)>a 2 (k)>...>d' Ns (k)} 9 where a x {k) is the 
largest singular value estimate and & Ns (k) is the smallest singular value estimate for 
subband k. When the singular value estimates for each diagonal matrix S(fc) are 
ordered, the eigenvectors (or columns) of the associated matrices U(£) and V(£) are 
also ordered correspondingly. After the ordering, & x (k) represents the singular value 
estimate for the best eigenmode for subband k 9 which is also often referred to as the 
"principal" eigenmode. 

[0082] A "wideband" eigenmode may be defined as the set of same-order eigenmodes 

of all subbands after the ordering. Thus, the w-th wideband eigenmode includes the /ra- 
th eigenmodes of all subbands. Each wideband eigenmode is associated with a 
respective set of eigenvectors for all of the subbands. The "principal" wideband 
eigenmode is the one associated with the largest singular value estimate in the matrix 

E(&) for each of the subbands. 
[0083] The user terminal can transmit a steered reference on the uplink. The uplink 

steered reference for the m-th wideband eigenmode may be expressed as: 

Xu P ^W = K m (A:)y ut>/n (/:) /? ^) , for k e K , Eq (34) 
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where v m m (fc) is the m-th column of the matrix V ut (fc) for the fc-th subband, with 

Y ut (*) = [y uU (*) 5U 2 (*) - i^^C*)] ; and 

/?(£) is the pilot symbol for the k-th subband. 
[0084] The received uplink steered reference at the access point may be expressed as: 

Lup^ m (k) = 5,(1)^.(4) + n up (*) , for k e K , Eq (35) 

= H,(i)K il (t)4 i(-l (*)p(« + n up (*) 

-a^(*)*.(*)p(*)+5,(« 

where u apm (/:) is the m-th column of the matrix U ap (£) for the k-th subband, with 

U ap W = [u apa (/:) u ap2 (*) ... u ap ^(fc)];and 

^" m (^) is the singular value estimate for the fc-th subband of the m-th wideband 
eigenmode. 

[0085] The access point can obtain an initial matched filter matrix M ap (£) , for ks K , 

based on the uplink steered reference, as described above. The access point may 
thereafter obtain an enhanced orthogonal matched filter matrix M ap (fc), for ke K , 

based on M ap (A:) and using any one of the orthogonalization techniques described 
above. 

[0086] Using QR factorization, the matrix M ap (fc) may be obtained as: 

MlpW = Q ap (^)Sap(*) >or Eq(36a) 
^(*) = fc(*)9^(*) = . Eq(36b) 

where Q ap (£) is a unitary matrix that is the ortho-normal basis for M ap (£) ; 
R ap (A:) is a diagonal matrix derived based on M ap (£) ; and 
lap(*)=]Lp(*) and 0; p (*) = Q ap (*). 

[0087] Using mean square error computation, the matrix M ap (A:) may be obtained as: 
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K^) = %l p (k)U Map (k)\^ ap (k) = %l p (k)vl(k) JorkeK, Eq(37) 

where 

C ap (*) = l ap (*)M ap (*) = U Map (*)S Map (*)Vl p (*) ; for k e K , Eq (38) 
2 ap (&) is the diagonal matrix whose elements are the 2-norm of the columns of 

Ml p (k) ; and 
U a >) = UMa P (*)YL p (*) . 

[0088] Using polar decomposition, the matrix M ap (£) may be obtained as: 

Map(*) = £ P (*)Z ap (^ JorkGK. Eq(39) 

where 

C ap (fc) = l ap (*)M ap (*) = Z ap (*)P ap (*) , for * e K ; and Eq (40) 

u! p (fc) = z ap w . 

[0089] The matrix M ap (fc) may be used by the access point for matched filtering of 

uplink data transmission from the user terminal, as described below. 
[0090] The spatial processing performed by the user terminal to transmit data on 

multiple eigenmodes on the uplink may be expressed as: 

x up (*)4 ul (^)Y ut (i)s up (i) ,for*GK, Eq(41) 

where s up (£) is the data vector and x up (k) is the transmit vector for the fc-th subband 

for the uplink. Uplink data transmission can occur on any number of wideband 
eigenmodes from 1 to Afc. 
[0091] The received uplink data transmission at the access point may be expressed as: 

r up (fc) = H up (A:)x up (A:) + n up (/:) 

= H up (*)K m (*)Y m (*)s up (*) + n up (*) , for k €= K , Eq (42) 

= y ap (/:)i(A:)s up (A:) + n up (/:) 

where r^(k) is the receive vector for the uplink data transmission for the k-th 
subband. 

[0092] The matched filtering by the access point may be expressed as: 
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§u P (*) = G ap (*)M ap (*)r up (*) 

= G ap (*)!!„ (*)f£ (U ap (*)|(*)s up (k) + n up (*)) , for k e K , Eq (43) 

=§ U p(*) + n U p(*) 

where !(/:) = diag (a u (k) cr 2 2 (k) ... & Nt „ r (k)) ; and 

G ap (*) = diag (l/oF, 2 ,^) l/«? 2 2 2 (*) ... l/^ 7 Wr (*)) . 

[0093] For the TDD MEMO system, the access point may also use the matrices U ap (&) , 

for k e K , for spatial processing for data transmission on the downlink to the user 
terminal. The spatial processing performed by the access point to transmit data on 
multiple eigenmodes on the downlink may be expressed as: 

x dD W = K ap (A:)U; p (fe)s dn (^) JorkeK, Eq (44) 

where s dn (*) is the data vector and x^ik) is the transmit vector for the k-th subband 
for the downlink. Downlink data transmission can similarly occur on any number of 
wideband eigenmodes from 1 to Ns- 
[0094] The received downlink data transmission at the user terminal may be expressed 

as: 

r dn (A:) = H dn (/:)x dn (A:) + n dn (A:) 

= H dn (^)K ap (t)y; p (i)s dn (i) + ri dD (t) JorksK, Eq (45) 

«Y m (*)i(*)s dn (fc) + n dn W 
where r dn (fc) is the receive vector for the downlink data transmission for the k-th 
subband. 

[0095] The matched filtering by the user terminal may be expressed as: 

idn(*) = G ot (*)M UI (k^ik) 

= G ut (k)± T (k)V M (*)(¥„, (VUVSto (k) + (k)) , for k e K , Eq (46) 
= s dn (/:) + n dD (A:) 

where M ul (k) = £ (£)V ul (k) is the matched filter for the user terminal; 
i(*) = diag (<7 U (*) & 22 (k) ... & NsJHm (*)) ; and 
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G m (*) = diag (l/<7, 2 ,(*) l/<7 2 2 2 (*) ... l/& 2 NsNs (k)) . 

The diagonal matrix is derived from the singular value decomposition shown in 
equation (32). 

[0096] Table 1 summarizes the spatial processing at the access point and user terminal 

for both data transmission and reception on multiple wideband eigenmodes. 
Table 1 





Downlink 


Uplink 


Access 
Point 


Transmit : 

x dn (A:) = K ap (fc)U; p (fc)s dn a) 


Receive : 

§u P (*) = G ap (/:)r(/:)U a W p (A:)r up (*) 


User 
Terminal 


Receive : 

id. (*) = G ul (k)t T (k)V ut (k) LdQ (*) 


Transmit : 

x U p(^) = K ut (/:)V ut (*)s up (/:) 



In Table 1, s(k) is the data vector, x(*0 is the transmit vector, r(k) is the receive 
vector, and s(k) is an estimate of the data vector s(k) , where all vectors are for 
subband k. The subscripts "dn" and "up" for these vectors denote downlink and uplink 
transmissions, respectively. 

[0097] It can be shown that the use of the matrices U ap (&), for keK y (with 

orthogonalized columns) for spatial processing for downlink data transmission can 

provide substantial improvement over the use of matrices U ap (/:) , for k e K , (with un- 

orthogonalized columns) obtained from the initial matched filter matrices M ap (£), for 
*e K . 

[0098] FIG. 4 shows transmission of steered reference and data on the downlink and 

uplink for an exemplary transmission scheme. The MIMO pilot is transmitted on the 
downlink by the access point in each TDD frame (block 412). The user terminal 
receives and processes the downlink MIMO pilot to obtain an estimate the downlink 

channel response H^fc), for ks K . The user terminal then estimates the uplink 
channel response as H CTip (/:) = H cda (k) and performs singular value decomposition of 
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H cup (k) to obtain the matrices E(&) and V u ,(fc), for k e K , as shown in equation (32) 
(block 414). 

[0099] The user terminal then transmits the uplink steered reference on the RACH or 

the RCH using the matrices V ut (fc), for k e K , as shown in equation (34) and FIG. 3, 

during system access (step 422). The columns of V ut (£) are also referred to as steering 
vectors when used for data transmission. The access point receives and processes the 
uplink steered reference on the RACH or the RCH to obtain the matrices %(k) and 

U ap (&), for keK, as described above (step 424). The columns of U ap (fc) are 

eigenvectors that may be used for both data reception as well as data transmission. The 
user terminal may thereafter transmit the uplink steered reference and data on the RCH 

using the matrices V ut (£), for k e K , as shown in equation (41) and FIG. 3 (step 432). 
The access point receives and processes the uplink steered reference on the RCH to 

update the matrices E(£) and U ap (fc), for ke K (step 434). The access point also 
performs matched filtering for the received uplink data transmission using the matrices 
l(jfc) and U ap (*) (also step 434). 
[00100] The access point may thereafter transmit an optional downlink steered reference 

and data on the FCH using the matrices U ap (fc) , for k e K , as shown in equation (44) 
and FIG. 3 (step 442). If a downlink steered reference is transmitted, then the user 
terminal can process the downlink steered reference to update the matrices %{k) and 
V ut (k) , for k s K (step 444) and may also perform orthogonalization to ensure that the 
columns of V m (fc) are orthogonal. The user terminal also performs matched filtering 

for the received downlink data transmission using the matrices %{k) and V ul (£) (also 
step 444). 

[00101] The pilot and data transmission scheme shown in FIG. 4 provides several 

advantages. First, the MIMO pilot transmitted by the access point may be used by 
multiple user terminals in the system to estimate the response of their respective MIMO 
channels. Second, the computation for the singular value decomposition of H (it), 
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for k e K , is distributed among the user terminals (i.e., each user terminal performs 
singular value decomposition of its own set of estimated channel response matrices for 

the usable subbands). Third, the access point can obtain the matrices %{k) and U ap (&) , 

for e K , which are used for uplink and downlink spatial processing, based on the 
steered reference without having to estimate the MIMO channel response. 

[00102] Various other transmission schemes may also be implemented for MIMO and 

MIMO-OFDM systems, and this is within the scope of the invention. For example, the 
MIMO pilot may be transmitted by the user terminal and the steered reference may be 
transmitted by the access point. 

[00103] FIG. 5 shows a block diagram of an embodiment of an access point 210x and a 

user terminal 220x in MIMO-OFDM system 200. For clarity, in this embodiment, 
access point 210x is equipped with four antennas that can be used for data transmission 
and reception, and user terminal 220x is also equipped with four antennas for data 
transmission/reception. In general, the access point and user terminal may each be 
equipped with any number of transmit antennas and any number of receive antennas. 

[00104] On the downlink, at access point 210x, a transmit (TX) data processor 514 

receives traffic data from a data source 512 and signaling and other data from a 
controller 530. TX data processor 514 formats, codes, interleaves, and modulates the 
data to provide modulation symbols, which are also referred to as data symbols. A TX 
spatial processor 520 then receives and multiplexes the data symbols with pilot 

symbols, performs the required spatial processing with the matrices U ap (&) , for k e K , 

and provides four streams of transmit symbols for the four transmit antennas. Each 
modulator (MOD) 522 receives and processes a respective transmit symbol stream to 
provide a corresponding downlink modulated signal. The four downlink modulated 
signals from modulators 522a through 522d are then transmitted from antennas 524a 
through 524d, respectively. 
[00105] At user terminal 220x, four antennas 552a through 552d receive the transmitted 

downlink modulated signals, and each antenna provides a received signal to a respective 
demodulator (DEMOD) 554. Each demodulator 554 performs processing 
complementary to that performed by modulator 522 and provides received symbols. A 
receive (RX) spatial processor 560 then performs matched filtering on the received 



030018 

EL977102456US 



29 



symbols from all demodulators 554a through 554d to provide recovered data symbols, 
which are estimates of the data symbols transmitted by the access point. An RX data 
processor 570 further processes (e.g., symbol demaps, deinterleaves, and decodes) the 
recovered data symbols to provide decoded data, which may be provided to a data sink 
572 for storage and/or a controller 580 for further processing. 
[00106] RX spatial processor 560 also processes the received pilot symbols to obtain an 

estimate of the downlink channel response, H cdn (fc), for £e K . Controller 580 may 
then decompose each matrix H cdn (fc) to obtain %{k) and V ut (fc). Controller 580 may 
further derive (1) the downlink matched filter matrices M ut (fc) , for /: e K , based on 
%{k) and V ut (A:), and (2) the scaling matrices G ut (£), for k e K , based on %{k) . 
Controller 580 may then provide M ut (&) to RX data processor 560 for downlink 

matched filtering and V ut (k) to a TX spatial processor 590. 

[00107] The processing for the uplink may be the same or different from the processing 

for the downlink. Data and signaling are processed (e.g., coded, interleaved, and 
modulated) by a TX data processor 588, multiplexed with pilot symbols, and further 
spatially processed by TX spatial processor 590 with the matrices V ut (£), for it e K . 
The transmit symbols from TX spatial processor 590 are further processed by 
modulators 554a through 554d to generate four uplink modulated signals, which are 
then transmitted via antennas 552a through 552d. 

[00108] At access point 510, the uplink modulated signals are received by antennas 524a 

through 524d and demodulated by demodulators 522a through 522d to provide received 
symbols for the uplink steered reference and data transmission. An RX spatial 
processor 540 then processes the received uplink steered reference to obtain estimates of 
n m CF m , for k e K and me {1 ... N s } , which are provided to controller 530. Controller 

then obtains M ap (fc) and Z(£) based on the estimates of u m cr m , performs 

orthogonalization of M ap (£) to obtain M ap (£) and U ap (fc), and derives G ap (fc) based 

on Z(k) . Controller 580 then provides M ap (fc) and G ap (£) to RX spatial processor 
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540 for uplink matched filtering and U ap (fc) to TX spatial processor 520 for downlink 
spatial processing. 

[00109] RX spatial processor 540 performs matched filtering of the received uplink data 

transmission with M ap (fc) and G ap (fc) to provide recovered data symbols, which are 

further processed by an RX data processor 542 to provide decoded data. The decoded 
data may be provided to a data sink 544 for storage and/or controller 530 for further 
processing. 

[00110] Controller 530 performs the processing to obtain the matched filter matrices 

M ap (fc) and the scaling matrices G ap (&), for ke K , for uplink data transmission and 

the matrices U ap (A:), for keK, for downlink data transmission. Controller 580 
performs the processing to obtain the matched filter matrices M ut (fc) and the scaling 
matrices G ut (fc), for ke K , for downlink data transmission and the matrices V ut (fc), 

for ke K , for uplink data transmission. Controllers 530 and 580 further control the 
operation of various processing units at the access point and user terminal, respectively. 
Memory units 532 and 582 store data and program codes used by controllers 530 and 
580, respectively. 

[00111] FIG. 6 shows a block diagram of the spatial processing performed by access 

point 210x and user terminal 220x to transmit data on multiple eigenmodes on the 
downlink and uplink. 

[00112] On the downlink, within TX spatial processor 520 at access point 210x, the data 

vector ^(fc) for each subband k is first multiplied with the matrix U ap (&) by a unit 
610 and further multiplied with the correction matrix K ap (£) by a unit 612 to obtain the 

transmit vector x<m(£) f° r subband k. The columns of the matrix U ap (fc) are 
orthogonalized as described above. The transmit vectors x^ik) , for ke K , are then 
processed by a transmit chain 614 within modulator 522 and transmitted over the 
MIMO channel to user terminal 220x. Unit 610 performs the spatial processing for 
downlink data transmission. 
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[00113] At user terminal 220x, the downlink modulated signals are processed by a 

receive chain 654 within demodulator 554 to obtain the receive vectors r^Ot), for 
ke K . Within RX spatial processor 560, the receive vector (k) for each subband k 

is first multiplied with the matched filter matrix M ut (fc) by a unit 656 and further 
multiplied with the scaling matrix G m (k) by a unit 658 to obtain the vector s^CA:), 
which is an estimate of the data vector (k) transmitted for subband k. Units 656 and 
658 perform the downlink matched filtering. 
[00114] On the uplink, within TX spatial processor 590 at user terminal 220x, the data 

vector s up (fc) for each subband k is first multiplied with the matrix V ut (fc) by a unit 

660 and then further multiplied with the correction matrix K ul (&) by a unit 662 to 
obtain the transmit vector x up (fc) for subband k. The transmit vectors x up (fc), for 

ke K , are then processed by a transmit chain 664 within modulator 554 and 
transmitted over the MIMO channel to access point 210x. Unit 660 performs the spatial 
processing for uplink data transmission. 
[00115] At access point 210x, the uplink modulated signals are processed by a receive 

chain 624 within demodulator 522 to obtain the receive vectors r^fc), for ke K . 

Within RX spatial processor 540, the receive vector r up (£) for each subband k is first 

multiplied with the matched filter matrix M ap (£) by a unit 626 and further multiplied 

by the scaling matrix G ap (fc) by a unit 628 to obtain the vector s up (fc), which is an 

estimate of the data vector s up (fc) transmitted for subband k. Units 626 and 628 

perform the uplink matched filtering. 
[00116] The techniques described herein to derive eigenvectors for spatial processing 

may be implemented by various means. For example, these techniques may be 
implemented in hardware, software, or a combination thereof. For a hardware 
implementation, the elements used for these techniques may be implemented within one 
or more application specific integrated circuits (ASICs), digital signal processors 
(DSPs), digital signal processing devices (DSPDs), programmable logic devices 
(PLDs), field programmable gate arrays (FPGAs), processors, controllers, micro- 
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controllers, microprocessors, other electronic units designed to perform the functions 
described herein, or a combination thereof. 
[00117] For a software implementation, the techniques may be implemented with 

modules (e.g., procedures, functions, and so on) that perform the functions described 
herein. The software codes may be stored in a memory unit (e.g., memory units 532 
and 582 in FIG. 5) and executed by a processor (e.g., controllers 530 and 580). The 
memory unit may be implemented within the processor or external to the processor, in 
which case it can be communicatively coupled to the processor via various means as is 
known in the art. 

[00118] Headings are included herein for reference and to aid in locating certain 

sections. These headings are not intended to limit the scope of the concepts described 
therein under, and these concepts may have applicability in other sections throughout 
the entire specification. 

[00119] The previous description of the disclosed embodiments is provided to enable any 

person skilled in the art to make or use the present invention. Various modifications to 
these embodiments will be readily apparent to those skilled in the art, and the generic 
principles defined herein may be applied to other embodiments without departing from 
the spirit or scope of the invention. Thus, the present invention is not intended to be 
limited to the embodiments shown herein but is to be accorded the widest scope 
consistent with the principles and novel features disclosed herein. 



WHAT IS CLAIMED IS: 



